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Theoretical limit of the minimal magnetization switching field and the optimal field 

pulse for Stoner particles 
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The theoretical limit of the minimal magnetization switching field and the optimal field pulse 
design for uniaxial Stoner particles are investigated. Two results are obtained. One is the existence 
of a theoretical limit of the smallest magnetic field out of all possible designs, ft is shown that 
the limit is proportional to the damping constant in the weak damping regime and approaches 
the Stoner- Wohlfarth (SW) limit at large damping. For a realistic damping constant, this limit is 
more than ten times smaller than that of so-called precessional magnetization reversal under a non- 
collinear static field. The other is on the optimal field pulse design: If the magnitude of a magnetic 
field does not change, but its direction can vary during a reversal process, there is an optimal design 
that gives the shortest switching time. The switching time depends on the field magnitude, damping 
constant, and magnetic anisotropy. However, the optimal pulse shape depends only on the damping 
constant. 

PACS numbers: 75.60.Jk, 75.75.+a, 85.70.Ay 
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Fabrication^], Q and manipulation^ of magnetic 
single-domain nano-particles (also called the Stoner par- 
ticles) are of great current interests in nano-technology 
and nano-sciences because of their importance in spin- 
tronics. Magnetization reversal, which is about how to 
switch a magnetization from one state to another, is an 
elementary operation. One important issue is how to 
switch a magnetization fast by using a small switching 
field. The switching field can be a laser light or a spin- 
polarized electric curre nt IH E| , or a magnetic field0,@. 



Many reversal schemes |9|, [10j have been proposed and 
examined. However, the issue of theoretical limits of the 
smallest switching field and the shortest switching time 
under all possible schemes are not known yet. Here we 
report two theorems on the magnetic-field induced mag- 
netization reversal for uniaxial Stoner particles. One is 
about the theoretical limit of the smallest possible switch- 
ing field. The other is about the optimal field pulse for 
the shortest switching time when the field magnitude is 
given. 

Magnetization M = rnM of a Stoner particle can be 
conveniently described by a polar angle 9 and an az- 
imuthal angle </>, shown in Fig.^, because its magnitude 
M does not change with time. The dynamics of magneti- 
zation unit direction in is governed by the dimensionlcss 
Landau-Lifshitz- Gilbert (LLG) equation[3j,|8|, 

(1 + a 2 )^^- — —to x h-t — am x (to x h t ), (1) 
at 

where a is a phenomenological damping constant whose 
typical value ranges from 0.01 to 0.22 for Co films 
The total field h t = h + hi comes from an applied field 
h and an internal field h\ = — V^w(to) due to the mag- 
netic anisotropic energy density w(rh). Different particle 
is characterized by different w(fh). In our analysis, we as- 
sume it uniaxial with the easy axis along the z-direction, 



u>(cos6*) and hi = — 



J* dw(cosO) 

~ d(cose) 



z = /(cos 9)z. 



According to Eq. (|TJ, each field generates two motions, 
a precession motion around the field and a damping mo- 
tion toward the field as shown in Fig. ^i. In terms of 
9 — (p, Eq. (JTJ can be rewritten as 3] , 

(1 + a 2 )9 — h t jj + ahg — af(cos9) sin 6?, 

(1 +a 2 )sin6»0 = ah^ - hg + /(cos0) sin 6*. (2) 

Here hg and h$ are the field components along eg— and 
e^— directions of to, respectively. 

The switching problem is as follows: In the absence 
of an external field, the particle has two stable states, 
Too (point A) and —Too (point B) along its easy axis as 
shown in Fig.^j. Initially, the magnetization is too, and 
the goal is to reverse it to —too by applying an external 
field. There are infinite number of paths that connect 
the initial and the target state. LI and L2 in Fig. ^ 
are two examples. Each of these paths can be used as 
a magnetization reversal route (path). For a given path, 
there are infinite number of ways of designing magnetic 
field pulses such that the magnetization will move along 
the path. Different design leads to different switching 
time. Let h L,s (t) be the magnetic field pulse of design s 
along magnetization reversal route L. To proceed, a few 
quantities must first be introduced. 

Definition of switching field H L s : The switching field 
H L ' S of design s along route L is defined to be the largest 
magnitude of h L > s (t) for all t, H L - S = max{\h L > s (t)\; Vt}. 

Definition of minimal switching field H L on reversal 
route L: The minimal switching field along route L is 
defined to be the smallest value of H L,S for all possible 
designs s that will force the magnetization to move along 
L, i.e. H L = min{H L ' s -ys}. 

Definition of theoretical limit of minimal switching field 
H c : The switching field limit H c is defined as the small- 
est value of H L out of all possible routes, i.e. H c = 
min{H L ;yL}. 
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FIG. 1: a, Two motions of magnetization m under field h: 
—in x h and — rh x (m x h) describe the precession and dissi- 
pation motions, respectively, b, Points A and B represent the 
initial and the target states, respectively. The solid curve LI 
and dashed curve L2 illustrate two possible reversal routes. 

If the applied field is restricted to be static, H L forms 
the so-called modified Stoner-Wohlfarth (SW) astroid0, 
and H c is below, but close to, the SW field when the 
applied field is about 135° from the easy-axis 7J. The 
first issue is about the existence of H c . 

Theorem 1: For a given uniaxial magnetic anisotropy 
specified by /(cos#) = — d Q^ose) ' ^ ne theoretical limit 
of the minimal switching field is given by H c = —p===Q, 
where Q = max{f (cos 9) sin 9}, 9 G [0, ir}. 

Proof: To find the lowest possible switching field, it 
should be noticed that field along the radius direction h r 
of an external field does not appear in Eq. J5J . Thus one 
can lower the switching field by always putting h r = 0, 

and the magnitude of the external field is h = \Jh 2 g + h 2 ^. 

According to Eq. (0) , 9 and 4> are fully determined by hg 
and htj, and vice versa. It can be shown that h 2 can be 
expressed in terms of 9, (j), 9, and <f> 

g = h 2 =(1 + a 2 )9 2 + 2af(cos9) sin (9(9 

+ {a sin 9 <j)) 2 + sin 2 9[<j)~ f (cos 9)} 2 . (3) 

Here g(9, 9, <j>) does not depend explicitly on <j) for a uni- 
axial model. 

In order to find the minimum of g, it can be shown 
that 4> must obey the following equation: 

<j ) = f(cos9)/(l + a 2 ), (4) 

which is from = and gf| ((M) > 0. 

Eq. Q is a necessary condition for the smallest min- 
imal switching field. This can be understood as follows. 
Assume H c is the minimal switching field along reversal 
path L described by 9(t) = 9\(t) and cj)(t) = (f>i(t) (i.e. 
H c is the maximum magnitude of the external field that 
generates the motion of 9\ (t) and (j>i (t))- If 4>i (t) does not 
satisfy Eq. Q}, then one can construct another reversal 
path L* specified by 9(t) = 9\(t) and 4>(t) = <j)2(t), where 
(f>2(t) satisfies Eq. (0}. Because 9(t) and 9 are exactly 



the same on both paths L and L* at an arbitrary time 
t, the values of g(t) shall be smaller on L* than those on 
L at any t. Thus, the maximum g* — (H*) 2 of g on L* 
will be also smaller than that (H 2 ) on L, i.e. H* < H c . 
But this is in contradiction with the assumption that H c 
is the theoretical limit of the minimal switching field. 
Hence, <j)(t) must obey Eq. (0J on the optimal path that 
generates the smallest switching field, H c . 
Substituting Eq. |@J into Eq. ©, we have 

h 2 = yi + a 2 9 + af(cos9) sin 0/y/l + a 2 } 2 . (5) 

In order to make h minimum, one notices that Eq. © 
does not contain 4>. Thus optimal reversal path should 
satisfy 9 > during its reversal from 9 = to 9 — 
ir. Furthermore, the second term in Eq. © is anti- 
symmetric about 9 = 7r/2 plane (positive in the upper 
hemisphere and negative in the lower hemisphere). Thus, 
the largest value of h on the optimal reversal path oc- 
curs in the upper hemisphere. It equals J^ a i ; where 
Q = max {f (cos 9) sin 9}. — QED 

To have a better picture about what this theoretical 
limit H c is, we consider a well-studied uniaxial model, 
w(fh) = —km 2 /2, or / = kcos9. It is easy to show that 
the largest h is at 9 = 7r/4 so that Q = fc/2, and 



At small damping, H c is proportional to the damping 
constant. The result in the limit of a — > coincides 
with the switching field in reference 9 where the time- 
dependent field always follows the motion of magneti- 
zation. At the large damping, H c approaches the SW 
field[j| when a non-collinear static switching field is 135° 
from the easy-axis. The solid curve in Fig. [21 is H c ver- 
sus a. For comparison, the minimal switching fields from 
other reversal schemes are also plotted. The dotted line 
is the minimal switching field when the applied field is al- 
ways parallel to the motion of the magnetization^. The 
curve in square symbols is the minimal switching field 
when a circularly polarized microwave at optimal fre- 
quencies is applied |9j. The dashed line is minimal switch- 
ing field under a non-collinear static field of 135° to the 
easy-axis. It saturates to the SW field beyond a c 0ll|- 

Although the theoretical limit of the switching field is 
academically important because it provides a low bound 
to the switching field so that one can use the theorem to 
evaluate the quality of one particular strategy, a design 
using a field at the theoretical limit would not be interest- 
ing from a practical point of view because the switching 
time would be infinite long. Thus, it is more important 
to design a reversal path and a field pulse such that the 
reversal time is the shortest when the field magnitude 
H (H > H c ) is given. An exact result is given by the 
following theorem. 
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Theorem 2: Suppose a field magnitude H does not 
depend on time and H > H c . The optimal reversal path 
(connects 9 = and 9 = ir) that gives the shortest switch- 
ing time is the magnetization trajectory generated by the 
following field pulse h(t), 

K{t) = 0, 

hg{t) =aH/Vl + a 2 , (7) 
h${t) = H/^/l + a 2 = he/a. 

Proof: The reversal time from A to B (Fig^ 3 ) is T = 
f£ dO/0. According to Eq. (J2J, one needs (h$ + ahe) to 
be as large as possible in order to make 9 maximal at 
an arbitrary 9. Since H 2 = h 2 + hg + hi, one has the 
following identity: 

(l + a 2 )H 2 = {l + a 2 )h 2 r + (h^ + ah e ) 2 + (h e -ah lj) ) 2 . (8) 

Thus, (hi), + ahg) reaches the maximum of vl + a 2 H 
when h r = and h$ = ah^, which lead to Eq. J7J), are 
satisfied. — QED 

Under the optimal design of 10, cf)(t) and 9(t) satisfy, 
respectively, Eq. 10} and 

9 = H/y/l + a 2 - a f (cos 9) sin (9/(1 + a 2 ). (9) 

For uniaxial magnetic anisotropy w(rh) = — km 2 /2, it 
is straight forward to integrate Eq. and to find the 
reversal time T from A to B (Fig^t) , 

T = j - ( " 2 + 1) " (10) 
k ^4(a 2 + l)H 2 /k 2 - a 2 

In the weak damping limit a — > 0, T m n/H while in the 
large damping limit a — > oo, T « , aTr — > oo. For 

the large field H — > oo, T « -J^+i-n ^ mverse ]y propor- 
tional to the field strength. Thus, it is better to make a 
as small as possible. Then the critical field is low, and 
the speed is fast (T ~ tt/H). Fig. EH shows the field de- 
pendence of the switching time for a = 0.1, where T and 
H are in the units of 2/k and k/2, respectively. 



How much could the so-called ballistic (precessional) 
reversal strategy 00 be improved? To answer the ques- 
tion, let us compare the switching field and time in the 
ballistic reversal with those of theoretical limits for uni- 
axial magnetic anisotropy w(rh) = —km 2 /2 and a = 0.1. 
According to reference 8, the smallest switching field (in 
unit of k/2) for the ballistic connection between A and B 
(Fig^j) is H = 1.02 applied in 97.7° to the easy z-axis, 
and the corresponding ballistic reversal time (in unit of 
2/k) is T = 5.87. On the other hand, the theoretical limit 
for the minimal switching field is H c w 0.1 from Eq. JHJ, 
about one tenth of the minimal switching field in the 
ballistic reversal 8J. For realistic value of a of order of 
0.01, the difference between experimentally achieved low 
switching field and the theoretical limit is of the order of 
hundred times. Thus there is a very large room for an im- 
provement. It is also possible to switch a magnetization 
faster than that of the conventional ballistic reversal by 
using a smaller field. For example, To achieve a reversal 
time of T = 5.87 along the optimal route, the field magni- 
tude can be as lower as H = 0.547 (instead of H = 1.02) 
according to Fig. El 

The field pulse given in Eq. Q requires a constant ad- 
justment of field-direction during the magnetization re- 
versal. To have a better ideal about the type of fields 
required, we plot in Fig. 0Ji the time dependence of x-, 
y- and z-component of the field while its magnitude is 
kept at H = 0.547. The time dependence of 9 and <f> is 
also plotted in Fig. ^d-c. 

Although the Stoner-Wohlfarth problem of magnetiza- 
tion reversal for a uniaxial model is of great relevance to 
the magnetic nano-particles, it is interesting to generalize 
the results to the non-uniaxial cases. So far, our results 
is on the magnetic field induced magnetization reversal, 
it will be extremely important to generalize the results to 
the spin-torque induced magnetization reversal. It should 
also be pointed out that it is an experimental challenge 
to create a time-dependent field pulse given by Eq. (Q) 
in order to implement the optimal design reported here. 

In conclusion, the theoretical limit of the magnetiza- 
tion switching field for uniaxial Stoner particles is ob- 



4 




FIG. 4: Time dependence of different field components, 9 and 
4> for uniaxial magnetic anisotropy w(m) = —kml/2 with 
a — 0.1 and H = 0.547 when the reversal path is optimal. 
The reversal time is T = 5.87. a, X-, y- and z-components of 
magnetic field, b, 9(t). c, <j>(t). 

tained. The limit is proportional to the damping con- 
stant at weak damping and approaches the SW field at 
large damping. When the field magnitude is kept to a 
constant, and the field direction is allow to vary, the op- 
timal field pulse and reversal time are obtained. 
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